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1. Introduction 

The robustness of the renowned result rj/s = l/Air is widely thought to be related 
to universal properties of black hole horizons. For p + 1-dimensional conformal field 
theories with a conserved charge that admit a dual holographic description a sim- 
ilar relation has been proposed for the diffusion constant, D, [1]. Recently, some 
understanding of how both results stated above could be derived from gravitational 
data that live at the black hole horizon has been pursued. In [2], by relating the 
charge diffusion in the boundary field theory to the analog process that occurs on the 
"stretched horizon" in the so called membrane paradigm, a closed analytic expression 
for the diffusion constant, D, was obtained as a a product of two factors. A local 
one, evaluated at the horizon, and a non-local one, involving an integral along the 
full range of the holographic radial direction. Certainly this smelled of holography, 
a point that was clarified later in [3], where the same formula was derived from a 
purely AdS/CFT construction. More recently [4] this factorized structure for the 
diffusion constant, was seen to match with the Einstein relation D = <tx~^- In it, 
the conductivity a was given by the local factor. The other one, namely the integral 
along r, corresponded exactly to the inverse of the charge susceptibility x~^y defined 
in equilibrium thermodynamics as the response of the charge density Ug to a change 
in chemical potential /i 
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A holographic model where we can compute the above quantities in a controlled 
manner is naturally given by a set of Nf "flavor" probe Dq branes, placed in the 
gravitational black hole background created by Nc 3> Nf "color" Dp branes. In order 
to model a chemical potential of the boundary field theory the world- volume U{1) 
gauge field At on the probe branes has to be switched on [6, 7]. These background 
ingredients are enough to compute the susceptibility (1.1). On the other hand, 
both the conductivity and the diffusion constant are transport coefficients, and their 
computation typically involves fluctuations. From the point of view of the underlying 
QFT, a and D are calculated from different channels of the retarded 2-point function 
of a conserved current, the transverse and longitudinal channels respectively. Up 
to Lorentz structures, the retarded current-current 2-point function is expressed in 
terms of two scalar functions Il±{k) and n||(fc). Pick's law for the conserved charge 
density implies the existence of a universal pole in the hydrodynamic limit k ^ 
in n||(fc) from where the diffusive dispersion relation can be read off 

w = -iDq"^ + 0{q^) . (1.2) 

On the other hand, from linear response theory we can calculate the conductivity 
from the following static limit ^ (see [1, 8]). 

CT = -Iim 1.3 

The fact that these three quantities obey the Einstein relation Dx = <y provides a 
non-trivial consistency check of the hydrodynamic picture of holography. Also the 
probe approximation is consistent with it. In fact, the parametric scaling Nf <^ Nc 
implicit in the probe approximation is correctly accounted for. Both x and a arise 
from normalized quantities that ultimately derive from the action, and therefore show 
the correct scaling NjN^ <^ N^ with respect to bulk quantities. On the contrary D 
stems from a pole condition and, at probe-level, does not scale with Nf. Backre- 
action should account for corrections of order Nf/Nc- This similar to the situation 
considered in [5] for the quotient ri/s 

The present paper is organized as follows. In section 2 we establish the class of 
models we shall be working with. In section 3 we will obtain an expression for the 
conductivity which includes and generalizes the one recently found in [9]. Next, in 
section 4, we study the charge susceptibility for which we also find an integral 
expression involving background fields that naturally extend the = case to 
non-zero baryon number. Finally, in 5, we turn to the diffusion constant, D. First 
we reanalyze the case of zero baryon number, following a simplified version of the 
argument in [3]. After that, we turn to the full system with chemical potential. The 
main difficulty stems from the fact that, in computing IIn, longitudinal fluctuations 
mix up with scalar perturbations of the probe brane profile. We prove that, unlike 
the case Uq = 0, the diffusion constant in general is not given by the one found in 

^For an electromagnetic conductivity, we should supplement this expression with a factor of 
the electromagnetic coupling. The same factor would show up in the susceptibility (see discussion 
in [1]) and they would cancel out in the diffusion constant. We will simply set e = 1.) 
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[2] using the membrane paradigm. So far we are able to provide a closed expression 
for D in terms of fluctuations, though not purely in terms of background quantities. 
Nevertheless, in section 6 we present numerical evidence that the Einstein relation 
holds also for rig 7^ 0. Assuming this is also true in general, the upshot is that the 
diffusion constant D can be computed from background fields as <JX~^- 



2. The holographic setup 

In this section we shall define the relevant notation. We will have in mind a generic 
holographic construction, where Nf Dq branes are embedded in the ambient metric 
of Nc ^ Nf Dp branes in the quenched approximation, where no backreaction of the 
flavor branes is taken into account. Therefore we shall parameterize a 10-dimensional 
metric in the form 

ds^ = gooir)dxl + 5'n(r)dXp + grrir)dr'^ + geeir)dnl + g^^{r)d'ilj'^ + g^^{r)dn'^_p_^ , 

(2.1) 

written to simplify the notation of a Dp/Dq intersection where the probe Dq branes 
wrap an n-sphere and the transverse space is spanned by coordinates ip, ip\ The 
presence of a horizon in the background metric aX r = th will be encoded in the 
following relations 

^00 = -gii{r)f{r) , 
grr = G{r)fir)~\ 

and /(r) = (r — ri{)F{r) with F{r) analytic at r = th- The embedding of the 
probe brane and the U{\) gauge field strength are parametrized by functions il){r) 
and A^j,{r) respectively, whose equations of motion are derived from the Born-Infeld 
action 

Sdbi = -NfTogJ rf^+i^e" V- det((7,5 + 27ia'Fab) • (2.2) 

The Wess-Zumino term plays no role in the following. Here we have introduced the 
tensor 

lab = gab + ^na'Fab , (2.3) 

and we shall be interested in solutions that involve a non-zero value for the temporal 
component of the gauge field 

7or- = 27ra'A[{r) = -7^0 , (2.4) 

corresponding to a non-zero chemical potential. For a nonvanishing value of the 
gauge potential as in (2.4) the matrix ■jab is non-diagonal, so we list here the inverse 
components (7"^7ca = ^c) 

7oo = -liif{r) ; 7rr = G(r)/(r)~^ ; /(r) = (r - rH)F{r) , 

7007rr- + (70r ) li: 
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where 7jj, ^qq, 7or and G{r) are regular at the horizon. The temperature can be 
expressed in terms of these coefficients in the near horizon hmit as 



rj. 1 liif 



4vr V-7oo7r 



(2.5) 



In general, these coefficients depend on the radial coordinate through the background 
proffie ip{r) and the world- volume gauge field Afj,{r). Reduced to such degrees of 
freedom, and integrated along the internal unit n-sphere of volume the DBI 
action acquires the form 



Sdbi = -NjTn,n^j (7oo7rr + ll) if.l^e , 

The charge density Ug can be obtained from the electric displacement [7] 



(2.6) 



with 7 = det7a;, without the angular components of the unit n-sphere and M = 
NfT£)^{27ia')'^Qn- From the equation of motion for At, it follows that the charge 
density is independent of the radial direction, drUq = 0. 

In the following sections we shall deal with the computation of retarded correlators. 
For this we consider fiuctuations in the worldvolume fields of the following form 

tlj{r,x) ipir) + ee-^('^^'°-9^'')^(r) , 
A,{r, x) ^ A,{r) + ee-^(-^'°-'?^'U^(r) . (2.8) 

If we expand the Dirac-Born-Infeld lagrangian to second order in powers of e 

C = Co + eCi + e^C2 + ... , (2.9) 

upon imposing the equations of motion for the background fields, the linear term, 
Ci, vanishes identically. The linearized equations for the perturbations and A are 
derived from the quadratic piece. In ref. [10] a detailed analysis of the retarded 
correlators was performed. We refer the reader to that reference for details on the 
equations of motion and the spectral functions. 



3. The conductivity 

The electrical DC conductivity may be obtained from the zero-frequency slope of the 
trace of the spectral function [8, 11] 



a 



lim 
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The spectral function is split into two orthogonal components (see [4, 10] for notation 
and conventions) 

X''^,{k) = -2{p - l)Imn^(fc) - 2Imnll(A;) . (3.2) 
Taking into account the boundary conditions for the scalar function Il-^{k) and nll(fc) 
n^(^,q = 0) = n||(a;,q = 0) ; n,, (^ = ±|q|) = . (3.3) 
From (3.1), the conductivity is uniquely obtained from the following expression 

Imn_|_(a;) , 
cr = - lim , 3.4 

independent of q (we shall confirm this independence shortly). In order to compute 
n_L we make use of the Minkowskian prescription given in [12] (see also [14]). To 
start with, we must solve the equations of motion for the transverse electric field 
fluctuations E± = ljA± with A± = ^2,3- 

El + dr log [e-^V^rY^] El - ^ E^ = . (3.5) 

The horizon r = rj^ is a regular singular point, and so we may perform a Frobe- 
nius analysis. Introducing as usual the dimensionless quantities vo = uj/{2tiT) and 
q = q/{2nT), the indicial exponents are found to be C± = ±«t»/2. We select C- 
corresponding to incoming wave boundary conditions and set 

E^= r^E^^reg. (3.6) 

In order to analyze the solution to (3.5) in the hydrodynamic limit, a power counting 
parameter A is introduced. Setting w — > Ato , q — >• Aq as well as 

- /-"" (<., + Ai5«„, + ...) = + A (<>„, - if <'„,log/) + ... 

(3.7) 

Expanding (3.5) in powers of A one easily solves for E'^^\,^g{r) in closed form 



dr 



The integrand diverg pole near the horizon. This enforces C2 = and we set 

C\ = \ which is simply an overall normalisation. This can then be used to solve to 
next order A in closed form 



<e,(-) = ,-,^^.^rr + ^ " (3-9) 



The logarithmic divergence at the horizon must be cancelled among the two contri- 
butions to E^^\r) by suitably tuning C4 = ivoC, where 



(3.10) 
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The finite part is fixed by adjusting C3 so that the boundary condition is E^±\eg{rH) = 
0. With this 

E^{r) = l + tw[ -^^^dr + OK,q2). (3.11) 

Once the exphcit solution has been found, the on shell boundary action can be 
evaluated 

AT " 



S 



B 



r=rB 



and from it, we obtain the correlator [12] 



-urrE'Ar)E^ir) 
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(3.12) 



(3.13) 



r^rg 



and we see from (3.11) that, to this order in A, there is no dependence on q. Therefore 
we can write the explicit expression for the conductivity as 



a 



1 , Imn^(t») 
lim 



2nT t»^o 



(3.14) 



In fact the imaginary part of (3.13) is independent of r, and therefore we may evaluate 
it at the horizon r — rather than at r^, because there we have explicit control 
over the singularities. Notice that E±{r) = 1 + 0{rx>), so finally we obtain 



a 



C 
2^ 



(3.15) 



where use of (2.5) and (3.10) has been made. For vanishing baryon number = 
the matrix 7^^ = Qah is diagonal and this expression fully agrees with that provided 
in [9]. It generalizes that result to the present context of holographic flavor systems 
with chemical potential. 



4. The susceptibility 

Let us now consider the susceptibility. This is an equilibrium quantity, given by the 
thermodynamical definition (1.1). On the one hand we have 



/i = At{rB) - Atiru) 



A[{r)dr 



(4.1) 



since A[rjj) = to have a well defined one-form at the horizon (see [7]). With this 

/ dA'(r) 

X={ I —^dr . (4.2) 



On the other, from (2.7) we can relate the charge density to A^{r) 

e~^^/HA[ 



Ug = J\f 



V- iloolrr + {2na'fA?) ' 



(4.3) 
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with H 



1\ 



. This equation can be inverted for A[ as a function of Ug 



a: 



riq 



-7007rr 



(4.4) 



Notice that in (4.2) we have emphasized the total derivative. This is so because 
on top of the exphcit dependence of A[{r) on n^, there is also a hidden one inside 
7rr = Qtt + i^' ijY gi>ip{'^) as well as in the factor e~'^^H = e~^'^7fj7gg,, stemming from 
the fact that the brane embedding ip{r) depends itself parametrically on Ug.^ In 
order to simplify (4.4) it is useful to make use of (2.7) and the following relations 



H 



1 



llr 



. (27ra 



/\2 



lOOlrr 

After some algebra one arrives at 

1 



(4.5) 



X 
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A 



dip' _ dtp 



dnn 



dn„ 



with 



(4.6) 



(4.7) 



It is remarkable to find the same combinations A and S that appear naturally in the 
equations of motion for the fluctuations (see later in (5.8 and ref. [10]). 



5. Charge diffusion 

In [2] the diffusion of a generic conserved charge was examined from the point of 
view of the membrane paradigm. A closed formula for this quantity was presented, 
and later reobtained in the context of the AdS/CFT correspondence in [3]. As a 
warm up exercise, we will rederive this expression in a somewhat simpler way here, 
first for Ug = 0. The expression obtained will also be valid for 7^ in the case of 
massless flavors. 



5.1 Zero baryon density 



The charge diffusion constant appears as a pole in the longitudinal correlator. In 
order to find the pole in it, one must study the equation of motion for the gauge 
invariant perturbations of the form E\\ = qAo + uAi. The relevant equation of 
motion is 



^11 ^11 







(5.1) 



^In [4] this fact was overlooked. Remarkably for vanishing baryon number it makes no difference 
since the factors multiplying dip/dug and dip' /duq vanish at Uq = 0. 



-7- 



From the diffusion pole we expect a dispersion relation of the form uj = —iDq^ + .... 
Therefore the natural hydrodynamic scaling is given in terms of a variable A as 
follows uj — s> A^cu, q — > \q. After checking that the indices near the horizon are again 
^-t = ±ito/2, a consistent expansion is 



E,{r) = f(rr-^ + X-Efl^ + ... = + A^^^ + ... (5.2 



,reg Hi^'efl 

Performing a Taylor expansion of the regular part around the horizon E\\^reg{f) = 

(0)' 

Ew^regi^H) + {r — ri^)£'|| j.^g{rH) + ••• we can obtain E^^ {th) solving the equation of 
motion (5.1) iteratively 

E'^\rH) = 1 , E^^m) = -~f\rH). (5.3) 

Moreover, inserting the ansatz (5.2) into the equation of motion (5.1) we may solve 
at lowest order in closed form 

Ef\r) = 1 - ,C?^ r — , (5.4) 



where the constant in front of the integral was fixed in accordance with (5.3) and 
(3.10). The relevant Green's function is proportional to Hn ~ limr_»rs E'^\lyr) j E\\{r b) 
where the boundary is at [4]. Hence the dispersion relation comes from demanding 
that E\\{rB) = 0. This provides us with the sought after dispersion relation w = 
—iD^q^ , with 

f'B Aj, 

^0 = e-*V^rV^7^^ / (5.5) 



Also here, for vanishing baryon number = ^ ^ab = 9ab diagonal and (5.5) 
coincides with the one in [2]. 

5.2 Finite baryon density 

In [13] a formula close to (5.5) (but not quite) was assumed to compute the value 
of Dq also at finite baryon density. The influence of Uq was supposed to show up 
through the dependence of the coefficients on the profile ip{r), which itself depends 
parametrically on the quark density. We shall see in this section that the answer that 
comes from examining the pole of the longitudinal propagator is more subtle. 
In the presence of a background value for At, the longitudinal perturbations exam- 
ined in the previous section mix with the scalar fiuctuations on top of the probe 
brane profile \E'(r). The coupled set of equations of motion can be written in the 
following form 

E[[ + ltiE[| + + €1$" + + = , (5.6) 

E[[ + %2E\\ + ^2^11 + €2*" + 5?2*' + «2* = , (5.7) 

where the explicit form of the coefficients is reproduced from reference [10] in ap- 
pendix A for completeness. From there one can easily see that all the coefficients 
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of \1/ are of order q. This implies that the natural variable is = would be 

natural to go with the gauge potentials A^). As before we define oj —>■ X^oj, q ^ Xq 
and it is clear that the coefficients of the terms multiplying \l/ are all of order A. 
Expanding to order A°, equation (5.6) becomes 

Or / 

+ log'(V=77"0'")<'' - (a*;„, + E*,„,) = . (6.8) 

Using (2.7) this equation can be integrated to give 

2TTa' 

By continuity in the limit — > we may reasonably expect that the constants 
Ci and C2 should be the same as in the previous subsection. A more rigorous 
derivation comes from comparing with the Frobenius expansion around, for instance, 
the horizon, as we did with (5.3) and (5.4). We have performed this comparison in 
the specific example of the D3/D7 and found perfect agreement. Therefore we finally 
set 

E}'\r) = 1+1 ^ A no,.. ■ (5.10) 



From the Dirichlet boundary condition E^^\rB) = we obtain the modification to 
the diffusion constant 

D = (5.11) 

2W /•(A^'(o) + 5^(o)) 



The sole purpose of presenting here this rather unwieldy formula for D is to empha- 
size that the influence of the baryon density on diffusion constant is not just given 
implicitely by modifying the embedding in Dq of (5.5). The coupling of the scalar 
modes that is only present at flnite density shows up as an explicit contribution^. It 
would be nice to have an expression for D in terms of purely background quantities. 
This would amount to an integration of the equation for the fluctuation \l/(r) which 
we have not achieved so far. However by providing numerical evidence that the Ein- 
stein relation holds (see section 6), we conjecture that indeed we can express D as 
a/x where both a and x ^i-re computable in terms of background quantities A^(r) 
and ^/^(r). 

Still an interesting case is that of massless quarks. For such a situation A and H 
vanish identically, and the diffusion constant is given hj D = Dq with the dependence 
on Uq coming through the dependence of 7afe on ^/'(r), and in this case the Einstein 
relation holds with flnite baryon density. 



■^A similar coupling can be seen to occur with the vector metric fluctutations in the case of the 
R-charged black hole. A parallel modification of the membrane-paradigm formula should correctly 
account for the diffusion at finite R-charge density, (see for example [15]). 
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6. Numerics, the D3/D7 case study 



In this section we shall focus on the case of a D3/D7 flavor setup. In this particular 
case we call the radial coordinate u, with uh = 1 and ub = 0. Now p = n = 3, = 
and 

7oo = -/ r ^ ; lu = - ; 7o« = -2WA n 6.1 

u u 

= ^ 4(1-V^>VH ' 7.. = i^(l-^M), (6.2) 



where 

/(«) = l-«2 ; T=^. (6.4) 

Here c? = nq{27ra')/M'rjj is a dimensionless parameter proportional to the density, 
whereas the coupling reads 

For finite baryon density, only results concerning the conductivity have been reliably 
established. In [16] a compact formula was found for the conductivity by imple- 
menting holographically Ohm's law on a probe D7 brane. A nontrivial check of this 
expression was provided in [10] by using the spectral functions and the same method 
explained here in section 3 



cr 



NcN/r 

An 



(1 - ^2)' + ■ (6.5) 



In figure 1 we show three plots for the three quantities, a, x and appropriately 
scaled, for different values of the baryon density d = 0.001, 0.005, 0.05, 0.2, 0.4 and 0.8. 
They are shown as a function of the dimensionless parameter m = M /T with M = 
2Mq/ V^, the mass gap. They all behave similarly, approaching different limiting 
values for large m. In the large m limit, a approaches an Uq dependent constant value 
whereas x dies off as ~ nq/m and D diverges proportionally to m after developing a 
minimum close to the d = curve. 

In the massless limit m — > the diffusion constant D reduces to Dq in (5.5) which 
can be expressed in terms of a hypergeometric function [20] 

1 3 4 



In figure 2 we have zoomed around the region of very small values of d. We observe 
the tri-valuedness that all of them exhibit for sufficiently low values oi d < d* = 
0.00315, which is related to the same property of the embeddings in this region of 
the parameter space [7, 13]. 
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Figure 1: In this figure we plot several curves for the diffusion constant D (in units 
of 1/2'kT), the susceptibility x (in units of NcNfT'^/2), and the conductivity a (in 
units of NcNfT/4ir), against the adimensional mass m = M /T. From bottom to top 
d = 0.001, 0.005, 0.05, 0.2, 0.4 and 0.8. 




The Einstein relation 

For the case of zero baryon number, the analysis performed in [4] confirmed the 
validity of the Einstein relation 

D=-. (6.7) 
X 

For nonvanishing baryon number this equation holds as well in the massless limit 
il){r) — s> 0. This can be seen from (3. 15), (4. 6), (5.5) and (5.11) after setting S = 
A = 0. For massive quarks, the validity of (6.7) has to be established numerically. 
D can be calculated both from a full numerical integration of (5.6) and (5.7) or by 
integrating ^(r) and using (5.11). We found complete agreement among them. Now, 
plotting D and (tx~^ leads to the set of curves shown in figure 3. We find a very 
good agreement and the discrepancies arise for high values of ipQ > 0.998 where the 
numerical computation of D is subject to large instabilities.^ 



^actually for D in figure 2, what we have plotted is the right hand side of (6.7). 
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Figure 3: Comparison of the diffusion constant D, and the quotient ajx- From bottom 
to top d = 0.005,0.05,0.2,0.4 and 0.8. We see remarkable agreement between the two 
calculations, up to high values of ipQ at which point the calculation of the excitations is 
subject to large numerical instabilities. 



7. Conclusions 

There is an intimate relationship between the membrane paradigm and the AdS/CFT 
prescription which is slowly unraveling and getting onto firmer grounds [3, 9, 18, 
17]. In this short note we followed the route of the first of these citations and 
worked fully within the AdS/CFT context. In this way we have generalized the 
closed formula obtained there for the diffusion constant D. We have also worked 
out the conductivity a and the susceptibility x- Fo^^ this last constant we provide 
an expression which matches the thermodynamic definition (1.1). An important 
technical detail in the calculation of x is tracking the implicit dependence of A'^ on 
Uq through the embedding profile ipijr) of the flavor brane. Including this contribution 
we have shown numerically that, at least for the D3/D7 case, the three constants 
obey the Einstein relation (6.7) also at finite baryon number 7^ 0. In the limit of 
massless quarks we have shown that this relation also holds in general. 
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A. Coefficients for coupled longitudinal system 



The coefficients that enter the system of differential equations (5.6) and (5.7) are 
reproduced here for completeness (see [10]). 
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(A.2) 
(A.3) 
(A.4) 
(A.5) 
(A.6) 



g(l-V''A) 



(A.7) 



(A.^ 



H(r) 



(A.9) 



with E and A as given in (4.7) and H (not that defined as in the main text) given 
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here by 

H(r) 



Yr{i - Tjj'^Y'g^i^) 



'^g^i>(y - i>'^Y''gi'i>) 



(A.IO) 
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